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A zigzag edge of a graphene nanoribbon supports localized zero modes, ignoring interactions. 
Based mainly on mean field arguments and numerical approaches, it has been suggested that inter- 
actions can produce a large magnetic moment on the edges. By considering the Hubbard model in 
the weak coupling limit, !/< t, for bearded as well as zigzag edges, we argue for such a magnetic 
state, based on Lieb's theorem. Projecting the Hubbard interactions onto the flat edge band, we 
then prove that resulting 1 dimensional model has a fully polarized ferromagnetic ground state. We 
also study excitons and the effects of second neighbor hopping as well as a potential energy term 
acting on the edge only, proposing a simple and possibly exact phase diagram with the magnetic 
moment varying smoothly to zero. Finally, we consider corrections of second order in U, arising 
from integrating out the gapless bulk Dirac excitations. 

One of the many fascinating properties predicted for graphene [1] is that a non-interacting nanoribbon with zigzag 
edges has bands of states with energy exponentially small in the ribbon width, localized at the edges [3] • Unzipping 
of carbon nanotubes has recently provided a technique for producing nano-ribbons with clean edges and Scanning 
Tunnelling Microscopy (STM) on such ribbons [5] has shown evidence for interacting edge states. It has been proposed, 
on the basis of mean field theory IH E] , density functional theory [7] and various numerical techniques [SHUT] , 
that Hubbard interactions may induce ferromagnetic order of the electrons in these edge states, with the moments on 
opposite edges ordering antiferromagnetically. Graphene edge magnetism looks promising for applications in nano- 
electronics |llj . However, there is no experimental evidence for edge magnetism in graphene ribbons and it is a 
matter of debate whether it is sufficiently robust to occur in realistic models [12) . thus motivating a deeper and 
more general understanding of its origins. One might try to regard this magnetism as an essentially one dimensional 
(ID) phenomenon, since it arises from edge states, but the ID Hubbard model is known to have a non-ferromagnetic 
ground state at all doping, being antiferromagnetic at half-filling. Recently p~3j [14], the effective ID model obtained by 
projecting the Hubbard interactions onto the edge states was studied numerically, and argued to lead to ferromagnetic 
order of an isolated edge. An interesting limit in which to try to prove edge magnetism is the weak interaction limit 
of the Hubbard model, U -C t, at half-filling. We take two steps towards proving edge magnetism in this limit. 
The first involves applying Lieb's theorem [15] to the contrasting cases of a nanoribbon with two zigzag edges (ZZ) 
versus a ribbon with one zigzag and one bearded edge (ZB). (See Figure 1.) Then we prove that the projected ID 
Hamiltonian has a fully polarized ferromagnetic ground state. We also study numerically other properties of the ID 
model, obtaining the electron or hole addition energy and showing that there are bound spin-1 excitons. We then 
consider two important particle-hole symmetry breaking perturbations: second neighbor hopping, ti (in the entire 
ribbon) and a potential energy, V e , acting on the edge atoms only, with the chemical potential maintained at the Dirac 
points of the bulk dispersion relation. Both perturbations lead |17j to the same new term in the ID Hamiltonian, 
oc t2 — V e = A. We argue that the fully polarized ground state survives up to a critical value of |A| of 0(U), beyond 
which the ground state may still be found exactly and has a smoothly decreasing edge magnetic moment. Integrating 
out the bulk excitations of the ribbon leads to both inter-edge and intra-edge interactions. The inter-edge interactions 
cx U 2 /(tW 2 ), where W is the ribbon width, produce antiferromagnetic order for the ZZ ribbon but ferromagnetic 
order for the ZB ribbon. Intra-edge interactions, cx U 2 /t, exhibit only a mild logarithmic singularity at low energies, 
arising from the gapless nature of the bulk Dirac excitations. 

We consider the Hubbard model at half-filling on a long ribbon of honeycomb lattice with periodic boundary 
conditions in the x-direction and zigzag edges. We actually find it convenient to first consider an upper zigzag edge 
and a lower bearded edge. [See Figure[l]] Let the number of atoms along the zigzag edge, of A type, be L. (Therefore 
the length of the ribbon is \/3aL where a is the nearest neighbour separation. We generally set -\/3a = 1.) L is also 
the number of hairs in the beard, ending at A sites. Noting that the number of A sites minus the number of B sites 
is L, it follows from Lieb's theorem [IS], that the ground states have spin S = L/2 for all values of U/t > 0. We 
label the width of the strip by another positive integer, W so that the width is 3aW/2 in the ZB case or a(3VF + 1) /2 
in the ZZ case. For the non-interacting model, U = 0, and W, L ^> 1, there are approximately L/3 zero energy 
states localized at the upper zigzag edge and 2L/3 localized at the lower bearded edge. The zigzag edge states have 
wave-vectors, 27r/3 < \k\ < n while the bearded edge states fill the rest of the Brillouin zone, |fc| < 27r/3. It follows 
from particle-hole symmetry that all of these states have exactly zero energy (for any W). These L zero modes are 
half- filled so we have a large degeneracy of ground states including a spin multiplet with spin L/2 obtained by filling 
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FIG. 1. A nanoribbon with an upper zigzag edge and lower bearded edge. In this example, L — 5 and W — 6. (We employ 
periodic boundary conditions in the ^-direction.) 



each of the L states with a single electron with the same polarization. Now consider turning on a very small positive 
U, in the case of large W and L. In the absence of boundaries we don't expect U to have a large effect on the ground 
state. A gapless non-magnetic Dirac liquid state is expected to persist up to a critical U c of 0(t). This follows from 
the fact that perturbation theory in U is infrared finite in the Dirac model; i.e. the 4-fermi interaction is irrelevant. 
Therefore, we expect the magnetic moment to live on the edges. The only physical explanation of the Lieb's theorem 
result seems to be that the fully polarized edge state multiplet persists as the (unique) ground state as we turn on U. 
This picture can be further substantiated by calculating the weak interaction between the upper and lower edges, of 
order U 2 / (tW 2 ) . This interaction is found to be ferromagnetic, as we show below. Thus there must be a spin sa L/6 
on the zigzag edge and « L/3 on the bearded edge, with these two spins coupled ferromagnetically. Now consider 
replacing the bearded edge by a second zigzag edge. Lieb's theorem |15j now implies a zero spin ground state since we 
have equal numbers of A and B sites. Now, at U = 0, we have approximately L/3 edge states on both lower and upper 
edges. These mix to form two bands, with 27r/3 < k < 47r/3, with energies exponentially small in W and symmetric 
around E = 0. Ignoring inter-edge interactions, we expect spin ss L/6 on both upper and lower edges. In this case, 
the intra-edge interaction of order U /(tW ) is antiferromagnetic, implying a zero spin ground state consistent with 
the result from Lieb's theorem. A further consistency check can be obtained by going smoothly between zigzag and 
bearded lower edges by turning on the hopping term, t' on the hairs. Lieb's theorem implies spin L/2 for all t, t' 
and U > 0. For t' — we have a ZZ ribbon together with L decoupled sites sitting below the lower edge. The ZZ 
ribbon has spin but we can obtain a state with total spin L/2 by polarizing the electron spins at the decoupled 
sites. Although the zigzag ribbon has total spin 0, for large W we expect that the upper and lower edges have spin 
ss L/6 with antiferromagnetic inter-edge coupling. Turning on t 1 produces an effectively antiferromagnetic coupling 
between the spin « L/6 on the lower zigzag edge and L/2 on the nearly decoupled sites. This gives a moment « L/3 
which is now ferromagnetically coupled to the upper edge, giving a total spin of L/2 as required by Lieb's theorem. 

Assuming that the ground state remains an unpolarized Dirac liquid up to U = U c , the magnetism of the edges 
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seems to follow from Lieb's theorem for large W. If the transition at U c is into a bulk antiferromagnetic state (with, 
for example, spin up on A sites and spin down on B sites) then the edge magnetism should persist, since it is of this 
type, and may be regarded as a sort of precursor of the bulk antiferromagnetic order. The simplest possibility is that 
the system goes from Dirac liquid into Mott-Hubbard insulator at U c but numerical evidence [TB] has been presented 
for a spin liquid phase at intermediate U, of unknown edge magnetic properties. 

We note that the above arguments also apply to carbon nanotubes [15) . Indeed, since we have been considering 
periodic boundary conditions in the ^-direction, we have actually been discussing tubes, of circumference L and 
length W. The magnetic moments exist on the upper and lower caps (i.e. rings) of the nanotubes with ferromagnetic 
or antiferromagnetic inter-ring coupling for a bearded or zigzag lower ring respectively. (The half-filled bulk of the 
nanotube might be in a ID version of a Mott-Hubbard insulating state but this only serves to weaken the effects of 
bulk states on edge states.) 

We may further substantiate this picture by considering [13) the weak intra-edge interactions of 0(U). Simply 
projecting the Hubbard interaction onto the zero energy states on the zigzag edge, in the large W limit, gives a 
Hamiltonian: 



n = \ E nk,k',q) 

k,k f ,q 



el(k + q)e a {k) - 5 q , Q 
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Here e a (k) annihilates an electron in an edge state with momentum k and spin a =f or J,. The interaction function 
is: 



T(k,k',q) = 
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The sum over k, k' and q is restricted to the band in which 27r/3 < k, k' , k + q, k' — q < 4-7r/3 and T(k,k',q) is 
strictly positive [16] . Periodic boundary conditions in the x-direction imply k = 2nn/L so the number of wave-vectors 
N ps L/3. Note that we are considering the case of half-filling in the entire lattice and that the edge Hamiltonian is 
therefore invariant under the particle-hole symmetry transformation: 



e a (k) <-» el(k). 



(3) 



This is highly unusual since normally a particle-hole symmetry transformation relates a particle and hole at different 
wave-vectors. Here with an exactly flat band, the particle and hole operators occur at the same wave-vectors. It is 
important to note that T(k,k',q) arose from summing the wave-function of the edge states over sites at arbitrary 
distance n from the zigzag edge and can be written: 



(4) 



where 



?n (fe) = e(l-|2cos||)Wl- (2 cos!) ( 2C ° S ^ 



(5) 



Thus, dropping the constant Eq, we may write: 



with 



n . q 



(6) 



On (q) = },gn[k)g n {k + q) 



q.O 



(7) 



It follows that all eigenstates of H are non-negative. It is can be seen that fully spin-polarized state is a zero energy 
eigenstate and therefore a ground state. To check this consider, for example, the representative fully polarized state 
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where all electron spins are in the up direction. Then clearly O n (q) annihilates this state for all non-zero q since the 
spin up terms in O n (q) try to produce a spin up electron in an occupied state with wave- vector k + q while the spin 
down terms try to annihilate a spin down electron in a vacant state of wave- vector k. O„(0) also annihilates this state 
since the occupancy of each single particle state is precisely 1. 

It is also possible, though more difficult, to argue that the fully polarized multiplet, of spin S = L/6, are the unique 
groundstates of the projected ID Hamiltonian. To prove that fully polarized states are the unique ground states of 
H we need to prove that the only states annihilated by O n (qYO n (q) for all n and q are fully polarized (that is, have 
maximal total spin). For convenience, in this paragraph, we take all momenta to be in the region of [— 7r/3, 7r/3], 
which can be obtained by shifting all of them by it. Suppose is such that for any n, q , O n {q)\tl>) = 0. Then we 
have 

Y^9n{k)g n {k + q) [el(k + q)e a (k) -<*,,„] \ip) =0 
fc 

= Y.9n(k) 9n (k + q) [ e ] T (k + q)e (7 (k) + ei(-k)e a (-k-q)-25 qfi ] \ip)=0 (8) 

fc>0 

since g n (k)g n (k + q) = g n (—k)g n (~k — q). (Repeated spin indices are summed in this section.) For fixed q and using 
the definition of g n we have, for any n that 

^(4sin(^)sin(^)y i |^)=0 (9) 



\^) = y/1 - (2 sin k/2) Vl - (2 sin(fc + g)/2) 2 [e\{k + q)e a {k) + e t (-fc)e CT (-A) - q) - 25 qfl ] \i>) (10) 

Since n runs from to oo, the number of independent momenta is L/3 and all the ^4 sin(|) sin(^i2)^ are different, 
the determinant of the Vandermonde matrix is non-zero so Eq. ([£]) is satisfied if and only if for any k, q we have 

[4(fc + q)e a (k) + et(-k)e a (-k - q) - 2S q , ] |V) - 0. (11) 



First, using Eq. (11) for q = 0, we get n(k) + n(—k) = 2, thus the only possible terms have n(k) = n(—k) = 1 
or n(k) = and n(—k) = 2 or vice versa. In general could be written as a linear combination of Fock states 
f[cj.(fc)|0 >. We first will show that in the expansion of |i/>) in terms of such states there is no Fock state which for 
any momentum k we have a vacancy or double occupancy in that momentum state. In other words, in the expansion 



of \ip), with condition (11 1 and n{k) +n(— k) = 2, only Fock states with singly occupied momentum states are allowed. 
Suppose that there is a state which has the property n(k) + n(—k) — 2 for any k and has double occupancy at 
momentum I (vacancy at momentum — I); call this state 4>: 

\4>) = \ — ,J±s — ,^tj—) (12) 
-i i 



and we suppose {(j)^) ^ 0. We impose the condition (11) on \tp) for q = —21 and k — I. Thus we should have 
ej.(— Z)e cr (/)|-0) = 0. Let us first look at the action of ej.(— l)e a (l) on <f), 

e t HK(0I4>> = 4HM0I ■ • • • • • > 4-t , • • •) 

i 

= 1 " ' -4-.' " ■ 'sj>' " ~ 1 " ■ 'sj>' " ' '4^' " ' } (13) 
-ii -ii 

Now, in order to satisfy the condition ej.(— Z)e CT (i)|^') = 0, we should have some other Fock states in the expansion 



of such that the action of ej.(— 2)e CT (Z) on them could cancel the terms created in the second line of Eq. (13) 
ej.(— l)e a (l) only acts on states with momentum I, —I thus does not change the spin configurations of the other (singly 
occupied) states. There are only three possible Fock states which have the same configurations of the singly occupied 
states: 

|i> = !■■■, 4 HK(Oli) = o 

-i i 
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|2> = | • • • ■ ■ ■ ■ • ■) 4 (-l)e a (l)\2) = | • • • . .. ..) 

-ll -i 1 

|3) = !•■■, ^.■■•^•••> 4HM0I3) = -|---,^,---,^,---> (i4) 



We see that none of these states are able to cancel the terms created on ( 13 ). Then the assumption is wrong and we 
have to have ((f>\ip) = 0. Having proven this, we show that in the expansion of \ip) in terms of singly occupied Fock 
states, only symmetric combinations like the following are acceptable 

I ■ • • , t . ' ' ' , I , ■ • •> + I ■ ■ • , I , ■ ■ ■ , t , ' ' •> (15) 

fcl k 2 fe'i fc 2 

Suppose that, the Fock expansion of \tp) has a term like | • • • , t > ' ' ' > I >'•")• Now by chosing k = ki and q = fc 2 — k\ 

fei fc 2 

we should have (eJ.(fc2)e CT (fci) + 4( — fci)e CT (— fca ) | V) — 0- We also have 

4(fc 2 K(fc!)l ■• • ■ •) = | • •• •• ••) 

fcl fc 2 fcl k 2 

eUh)e a (h)\ I t ,■■•> = -!••• • • • at , • • •) (16) 

fei fe 2 fc i fe 2 

Thus we see that the symmetric combinations leads to zero, while the anti-symmetric ones give us a non-zero result. 
Thus \ip) must be fully symmetric under exchanging all spins and is therefore of maximal spin. 

One interesting quantity that follows from the Hamiltonian is the energy to add a spin down electron of momentum 
k, which is the same as the energy to remove a spin up electron of momentum k: 

k' 

This quantity is plotted (at L — > oo) in Fig. [2j It vanishes linearly at the Dirac points, k — 2tt/3, 47r/3. In principle, 
e,t could be measured in Angular Resolved Photo-emission Spectroscopy (ARPES) experiments. The corresponding 
electron addition or removal energy is given by e/c. The corresponding density of states cx l/\dek/dk\ could be 
measured by Scanning Tunnelling Microscopy (STM). With a spin-polarized STM tip and an edge fully polarized in 
the z-direction, it would only be possible to tunnel in a spin-down electron or tunnel out a spin-up electron. 

We have calculated numerically the lowest energy particle- hole state of total momentum q, for L up to 602 (N = 200). 
This is plotted in Fig. [3] along with the bottom of the particle-hole continuum. We see that a strongly bound exciton 
exists for most values of q, as might be expected in this strongly interacting system. However, the binding energy 
vanishing at q = ±2it/3. This vanishing can be understood from the fact that £/, q) vanishes when k or fc' is at a 
band edge 27r/3 or 47r/3 so the zero energy particle and hole become non-interacting at wave-vector 27r/3 and —2tt/3 
or vice versa. 

While Lieb's theorem continues to imply a fully polarized ground state at sufficiently small U for any hopping terms 
between opposite sub-lattices (A to B), adding a small [0(£/)] second neighbor hopping term, t2, may destroy the 
fully polarized state. Likewise, a single site potential, V e , acting at the edge of the ribbon only, could destroy the fully 
polarized state. Temporarily ignoring interactions, the zigzag edge states survive at finite i 2 and V but develop a 
non-zero dispersion given, to first order in A = t^ — V e , by [T7]: e 2 (fc) —ep = A(2cos/c + l). breaking the particle-hole 
symmetry. Here we are assuming, for simplicity, that the bulk chemical potential is at the energy of the bulk Dirac 
points, which becomes ep = 3t 2 - (Shifting ep away from the Dirac points, the Hubbard interactions have a larger 
effect in the bulk rendering the edge model approach more questionable.) Including a small U, the energy to add a 
spin down electron or remove a spin up electron at momentum k now becomes 

E p/h (k) = e k ± A(2cosA; + l). (18) 

respectively, where tk is given in Eq. (17). E p (k) is plotted in Fig. [2] for several values of A. We see that for 
|A| < A c w 0.109J7, the energy to add an electron or hole remains positive, so the edge states remain undoped. A 
local minimum at k — tt develops in E p (k) for A > .087J7, and E p (k) becomes negative in the vicinity of k = n 
for A > A c . The lowest energy of a particle-hole state, and the bottom of the particle-hole continuum for various 
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FIG. 2. Energy to add a spin down electron of momentum k for various values of r\ = A/17 = (t2 — V e )/U. 




FIG. 3. Lowest energy particle-hole state (circles) and bottom of the particle-hole continuum (lines) for various values of 
77 = A/U = (t 2 - V e )/U. 
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FIG. 4. Edge magnetization versus A = ti — V e . 



values of A are shown in Fig. [3j We see that the exciton becomes unbound except for wave-vectors near zero, as 
|A| increases. For |A| > A c the edge states become doped, adsorbing electrons or holes from the bulk. The simplest 
assumption for |A| > A c is that a Fermi sea of spin-down electrons or spin up holes forms near k — tt, for A > A c or 
A < — A c respectively. This assumption is reasonable since there appear to be no bound excitons for A > A c . We 
also calculated for A near A c and L < 74, the lowest energy state with M = N/2 — 2, finding no states below the 
2-particle, 2-hole continuum, consistent with this assumption. These are exact eigenstates of the Hamiltonian of Eq. 
0, ( fl8| ) with the added holes or particles non-interacting. The corresponding exact result for magnetization versus 
A is plotted in Fig. [4j given this assumption. The non-interacting nature is a simple consequence of the fact that 
the on-site Hubbard model only gives interaction between electrons of opposite spin. On the other hand, we cannot 
rule out the possibility that the ground state for A > A c contains a finite density of spin-up holes as well as the 
spin-down electrons (and similarly for A < — A c ). In that case, Hubbard interactions have a non-trivial effect. In any 
event, adding a nearest neighbour Coulomb repulsion term to the bulk Hamiltonian has no effect on the projected 
edge Hamiltonian, since such a term acts between A and B sites whereas the zigzag edge states live entirely on one 
sublattice. On the other hand, a second neighbour Coulomb repulsion, U2, produces interactions between electrons 
with parallel spins in the projected edge Hamiltonian. While this doesn't change our conclusions qualitatively in the 
undoped phase, it will produce interaction effects in the ground state for the doped case even if contains only particles 
or only holes. However, we might expect U% <C in which case these effects could be quite small. Thus in general 
we expect a one or two component Luttinger liquid for |A| > A c . On the other hand the edge phase occurring for 
|A| < A c is definitely not a Luttinger liquid. Instead, it might be described as a fully spin-polarized semi- metal since 
all levels are filled with spin-up electrons and there is a non-zero electron and hole addition energy for all wave-vectors 
accept the band-edges, 2tt/3 and 4ir/3. 

There are also important effects of 0(U 2 /t) which arise from the interactions between bulk and edge states. We can 
consider integrating out the bulk states to obtain a low energy effective action for the edge states. Due to the gapless 
nature of the bulk Dirac spectrum, this produces long range retarded interactions among the edge excitations. Decay 
processes of edge into bulk electrons are forbidden by energy-momentum conservation but the Feynman diagrams of 
Fig. [5] induce quartic interaction terms. For large W and low energies we may calculate these interactions keeping 
only the low energy bulk states near the Dirac points, using the corresponding Dirac propagators. Note that we ignore 
interaction effects in the bulk, as discussed above. This is rather similar to an RKKY interaction. The interaction 



where 



involving the dynamical spin operators |20| Su / l(u> , q) on the upper and lower edge, respectively is: 

/dqduj 
-^-^Su(oj,q) ■ S L (-uj,-q)J intei (uj,q,W) (19) 

J inter ( W , q, W) = 2 U 2 J %^G{J, k, 0, W)G(u -<J,q- k, 0, W). (20) 

Here G(u, fc, 0, W) is the bulk free electron Green's function with momentum k in the x-direction at y = and y = W 
with appropriate zigzag or bearded boundary conditions and projected onto the sublattices corresponding to the upper 
and lower edge (A — A for zigzag-bearded or A — B for zigzag-zigzag) . Using the linearized, Dirac dispersion relation, 
which is valid at small 1/W, u/t and q, 

o,y . W) - & £ ^.^tfeAM (21) 

where e(fc) = wfI^I is the Dirac dispersion relation. Gzz is given by the same expression with u> replaced by ie(k x , k n ) 
in the numerator inside the sum. The sum over n can be taken up to an arbitrary ultra-violet cut-off whose value 
doesn't affect the behavior at small 1/W, uj/t and q. k n = im/W for the ZB case. Although the wave- vectors of edge 
modes are phase-shifted from these values in the ZZ case, this can be ignored at leading order in 1/W, allowing us to 
again use k n = irn/W. 

It is straightforward to evaluate J[ n i eT (u, q, W) numerically with the two types of edges. The characteristic scales 
for the uj and q dependence of </i n t e r are se ^ by t/W and 1/W, respectively. Since the energy scale of the inter- 
edge interaction in Eq. (21 1 is U 2 /(tW 2 ), it should be permissible to ignore the retardation, and evaluate Sj n ter a * 
lo = q = to calculate the properties of low energy states. This gives: 

JzB/zzm^Tc^^ (22) 
where the positive constant c is given by the convergent sums and integral: 

xfa ^ f°° n 2 m 2 1 

c =lA x Y (_!)»+". / dn—, ^ - « 0.20 (23) 

(A similar result was obtained in [J] for the ZZ case.) We see that the ground state for the zigzag-bearded ribbon 
has spin L/2 while that for the zigzag-zigzag case has spin 0, as shown above rigorously using Lieb's theorem. The 
remarkable fact that the change in sign of this tiny coupling drastically changes the spin of the ground state provides 
evidence for the polarized nature of the edge spins. There is also a large manifold of low energy states, which are 
simply the eigenstates of JzzSt • Sb with St — Sb — L/6 (in the ZZ case). 

Another important effect of 0(U 2 /t) is the intra-edge interaction, independent of W. For a zigzag edge by inte- 
grating out the low-energy bulk excitations, the spin part is: 



dqduj 

W) 



'intra = / 7^32 S ( w > «) ' 9Mntra( w > ?) ( 24 ) 



with 



J intra (uj,q) = 2U 2 [ d ^G(u;',k,y = y' = 0)G(u-u',q-k,y = y' = 0) (25) 

Jk,u>'/vf<A \ Z7r ) 

Now the free bulk Green's function, with zigzag edge boundary conditions, may be evaluated for a semi-infinite system, 
giving, at small k x (measured from a Dirac point) and small u>: 

<*>.*~v-*-Q~»ifg jSv^rkw (26) 



By using this green function, the ^ m |; ra of Eq. (25 1 is ultraviolet divergence and the integral should be cut off at some 
point A. Although the resulting t/ m ^ ra (w, q) is cut off dependent, by ignoring the weak retardation, the corrections to 
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FIG. 5. Feynman diagrams inducing edge interactions from integrating out bulk states. 
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the energy of the excitons is proportional to Ji n t Ta (0, l) ~ ^intra(^' wmcn is cut °ff independent. For small ui and 
q these 0(U 2 ) intra-edge interactions become more singular than the 0(U) terms by logarithmic factors of q 2 \n.q : 
giving for example the final equation in the main paper. 

As mentioned above, Eq. ( 25 ) only includes the effect of low energy bulk excitations; it is still possible that the high 
energy bulk excitations could wipe out this singularity of the exciton dispersion relation. By using the exact form of 
the bulk wavefunctions one can determine the exact intra-edge interaction of 0{U 2 ). This has a more complicated 
form than Eq. (25). Nonetheless, it can be shown that the only part of this interaction which contributes to this lng 
singular correction to the excition dispersion relation is the contribution from low energy bulk excitations of the form 
ofEq. ([25]). 

We leave a more detailed study of these effects of 0(U 2 ) and higher for the future. A reasonable approach might 
be to ignore the bulk interactions, since they are irrelevant, but analyse the bulk-edge Hubbard interactions using 
the renormalization group. This corresponds to a novel type of boundary critical phenomena in which the bulk is a 
massless (2+1) dimensional Dirac liquid and the edge is a one-dimensional spin-polarized semi-metal. The arguments 
based on Lieb's theorem imply that the edge magnetic moment remains stable against weak interactions. 

We would like to thank Ion Garate for his collaboration in the early stages of this project and Pawel Hawryak for 
interesting discussions. The research was supported in part by NSERC and ClfAR. We also thank the Galileo Galilei 
Institute where this work was completed. 
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